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I. INTRODUCTION 



Quantum states of a relativistic massive fermion are described by four-components wave func- 
tions called Dirac spinors. These wave functions, which are solutions of the Dirac equation, describe 
states of positive and negative energy. For the case where fermions carry the electric charge q, the 
electromagnetic interaction of fermions can be included by means of the electromagnetic fourpo- 
tential A 11 , which is introduced in the Dirac equation through the so called minimal substitution, 



changing the fourmomentum such as p^ 



qA^. It is also possible to introduce a linear 



harmonic potential in the Dirac equation by substituting p '— > p '— imufir, where m is the fermion 
mass, oj represents an oscillator frequency, r is the distance of the fermion respects to the origin of 
the linear potential and j3 = 7 corresponds to the diagonal Dirac matrix. 

The Dirac equation including the linear harmonic potential was initially studied by Ito et al. 

Cook and Ui et al. Q|. This system was latterly called by Moshinsky and Szczepaniak as 
Dirac oscillator [4J, because it behaves as an harmonic oscillator with a strong spin-orbit coupling 
in the non-relativistic limit. As a relativistic quantum mechanical system, the Dirac oscillator has 
been widely studied. Several properties from this system have been considered in (1+1), (2+1), 
(3+1) dimensions la]- 28]. Specifically, for the Dirac oscillator have been studied several properties 
as its covariance [6j], its energy spectrum, its corresponding eigenfunctions and the form of the 
electromagnetic potential associated with its interaction in (3+1) dimensions [7], its Lie Algebra 



symmetries 



the conditions for the existence of bound states 



its connection with supersym- 



metric (non-relativistic) quantum mechanics [10( , the absence of the Klein paradox in this system 



111 ], its conformal invariance its complete energy spectrum and its corresponding eigenfunc- 
tions in (2+1) dimensions [ljj, the existence of a physical picture for its interaction 14 1. For this 
system, other aspects have been also studied as the completeness of its eigenfunctions in (1+1) and 



n 



(3+1) dimensions |l5|], its thermodynamic properties in (1+1) dimensions [if]], the characteristics 



of its two-point Green functions 13], its energy spectrum in the presence of the Aharonov-Bohm 
effect 18j, the momenta representation of its exact solutions [19], the Lorenz deformed covari- 



ant algebra for the Dirac oscillator in (1+1) dimension [20f | . the properties of its propagator in 



(1+1) dimensions using the 
a Jaynes-Cummings model 



supersymmetric path integral formalism 



21[, its exact mapping onto 



221 ]. its nonrelativistic limit in (2+1) dimensions interpreted in terms 



of a Ramsey-interferometry effect 



23f ]. the existence of a chira. 



in (2+1) dimensions in presence of a constant magnetic field 



solutions using the Clifford algebra 



phase transition for this system 
241 ] . a new representation for its 



25[, its dynamics in presence of a two-component external field 
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261 ]. the relativistic Landau levels for this system in presence of a external magnetic field in (2+1) 



dimensions 



271 ] and its relationship with (Anti)-Jaynes-Cummings models in a (2+1) dimensional 



noncommutative space 



28]. 



Some possible applications of the Dirac oscillator have been developed. For instance, the 
hadronic spectrum has been studied using the two-body Dirac oscillator in [j], [3] and the ref- 
erences therein. The Dirac oscillator in (2+1) dimension has been used as a framework to study 
some condensed matter physical phenomena such as the study of electrons in two dimensional 
materials, which can be applied to study some aspects of the physics of graphene 30f|. This system 



has also been used in quantum optics to describe the interaction of atoms with electromagnetic 
fields in cavities (the Jaynes-Cummings model) 3o| . 

The standard point of view of the Quantum Field Theory (QFT) establishes that an excitation 
of one of the infinite degrees of freedom that constitute the free Dirac quantum field can be 
interpreted as a free relativistic massive fermion [31[. The free Dirac quantum field which describes 



the quantum dynamics of a non-interacting relativistic massive fermion can be seen as constituted 
by infinite decoupled quantum harmonic oscillators. In this scheme a non-interacting relativistic 
massive fermion is described as an excitation of a degree of freedom of the free Dirac quantum field 



311 ] . For the free Dirac field the energy quanta of the infinite harmonic oscillators are relativistic 



energies of free fermions or antifermions. An analogous situation is presented in the description 
of a non-interacting relativistic boson as an excitation of one of the infinite degrees of freedom 



(harmonic oscillators) that constitute a free bosonic quantum field 3lf] . For the free Boson field 



the energy quanta of the infinite harmonic oscillators are relativistic energies of free bosons. 

On the other hand, the Dirac oscillator, which describes the interacting system constituted by 
a relativistic massive fermion under the action of a linear harmonic potential, has not been studied 
as a quantum field up to now. In this direction, we show that it is possible to study consistently 
the Dirac oscillator as a QFT system. The main goal of this work is to consider the Dirac oscillator 
as a field and perform the canonical quantization procedure of the Dirac oscillator field in (1+1) 
and (3+1) dimensions. We can do it, in the first place, because the Dirac oscillator is characterized 

n 

by the absence of the Klein paradox [111 ] . that allows to distinguish between positive and negative 



energy states 



321 ]. In the second place, the Dirac oscillator is an interacting system in which there 



exist bound states that are defined by means of specific quantum numbers. This aspect means 
that the solution of the Dirac oscillator equation leads to both a well-known energy spectrum and 
an eigenfunction set which satisfies the completeness and orthonormality conditions in (1+1) and 
(3+1) dimensions [if]]. Thus the Dirac oscillator field can be consistently treated as a quantum 
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field because it can be written in terms of a Fourier expansion of that eigenfunction set. We observe 
that after the canonical quantization procedure, the Dirac oscillator field can be seen as constituted 
by infinite relativistic decoupled quantum linear harmonic oscillators. 

The canonical quantization procedure for the Dirac oscillator field starts from writing the Hamil- 
tonian operator for this system in terms of annihilation and creation operators satisfying the usual 
anticommutation relations. This procedure allow us to obtain the Feynman propagator for the 
Dirac oscillator field. Thus we are able to study the differences between the Dirac oscillator field 
and the free Dirac field. We note that while for the free Dirac field the energy quanta of the infinite 
harmonic oscillators are relativistic energies of free particles, for the Dirac oscillator quantum field 
the quanta are energies of relativistic linear harmonic oscillators. 

This work is presented as follows: First, in section 2 we introduce the notation used in this paper 
presenting the relevant aspects of the standard canonical quantization procedure for the free Dirac 
field in 3+1 dimensions. Next, in section 3 we study the Dirac oscillator system in 1+1 dimensions 
from a point of view of quantum relativistic mechanics; additionally we describe some aspects 
of the Dirac's sea picture for this system and we develop the canonical quantization procedure 
for the Dirac oscillator field in 1+1 dimensions. Posteriorly, in section 4 we study the canonical 
quantization procedure for the Dirac oscillator field in 3+1 dimensions. Finally, in section 5 we 
present some conclusions of this work. 



II. CANONICAL QUANTIZATION FOR THE FREE DIRAC FIELD 

In this section we present the main steps of the standard canonical quantization procedure of 



the free Dirac field following the procedure developed in 



311 ] . The equation of motion for a free 



relativistic fermion of mass m is given by the Dirac equation 

(n^-m)V(r,i) = 0, (1) 

where we have taken H = c = 1. In this equation 7^ represents the Dirac matrices obeying the 
Clifford algebra {7 /i ,7 l/ } = and ^j(f,t) represents the four-component spinor wave functions. 
The four linearly independent solutions from the Dirac equation are plane waves of the form 
if;l(f,t) = (2vr)- 3 / 2 Jj^Wripje-^^-P-^ , where E f is given by Ep = +^p 2 + m 2 . The two 
solutions of relativistic free-particle positive energy are described by r = 1, 2 while the two solutions 
of relativistic free-particle negative energy are described by r = 3, 4. The sign function e r takes 
the value e r = 1, for r = 1, 2, and e r = —1, for r = 3, 4. Additionally, the spinors w r {p) obey the 
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equation (7 M p M — e r m)w r (p) = 0. 

In a quantum field theory treatment, ip(r,t) is the free Dirac field. The Hamiltonian associated 
to this field is H = J d?x^ (—ia ■ V + m{3)ip, where the matrices a,f3 are defined as ai = 7°7*, 
/3 = 7°. The canonical quantization procedure of the free Dirac field starts by replacing the fields 
i/j(r, t) and ^(f, t) with the fields operators i/j(r, t) and tf(r, t) which obey the usual commutation 
relations of the Jordan- Wigner type [311 ] . The Fourier expansion of the free Dirac field operator 



ip(r,t), in terms of the plane waves functions, is written as 

where p-x = p^x^ = Egt — p- r and a(p, r) is an operator. The operators a(p, r) and its conjugated 
d^(p,r) satisfy the usual anticommutation relations. The Hamiltonian operator H associated to 
this system is defined in terms of the Dirac field operators ip(r, t) and t). Using the Fourier ex- 
pansion of the free Dirac field operator, H can be written as H = J d?pEp Ylr=i ™p,r + Sr=3 ™p,r > 
where r)a(p, r) is the particle-number operator and npf = a(p,r)a^ (p,r) is the 

antiparticle-number operator. In this Hamiltonian operator, which is positively defined, the Dirac's 
see picture has been used and the zero-point energy contribution has been subtracted 31]. 



Now it is possible to introduced the following canonical transformations over the annihilation 
and creation operators that allow to differentiate between operators associated to particles and 
antiparticles: c(p,+s) = a(p,l), c(p,—s) = a(p,2), <v(p,+s) = a(p,3), d^(p,—s) = a(p,A). With 
these transformations, the following identifications are possible: c(p,s) and &{p,s) represent the 
annihilation and the creation operators of particles, respectively; d(p,s) and d)(p,s) represent the 
annihilation and the creation operators of antiparticles, respectively. Using the operators c(p, s), 
c'(p, s), d(p, s) and d'(p, s), the Hamiltonian operator of the free Dirac field can be written as 31] 

s J 

where hp s = (p, s)c(p, s) represents the fermion-number operator and hp s = d) (p, s)d(p, s) rep- 
resents the antifermion-number operator. It is possible to observe that c'(p,s) creates a fermion 
of mass m having energy Eg, momentum p, charge +e and projection of spin +s, whereas cft(p, s) 
creates an antifermion of the same mass having the same energy, the same momentum, charge — e 
and projection of spin —s. In this scheme, a free relativistic massive fermion is described as an 
excitation from a degree of freedom of the free Dirac quantum field. For the free Dirac field, the 
energy quanta of the infinite harmonic oscillators are the relativistic energies Eg of free fermions 
or antifermions. 
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Finally, the Feynman propagator S^ /3 (x — y) of the free Dirac field is defined as iS^ /3 (x — y) = 
(0| T (^ a (x)^p(y)\ |0). By using the plane wave expansion, this propagator can be re-written as 
iS^g(x — y) = {ij^d^ — m)5 a /^iAF(x — y), where iAp(x — y) represents the Feynman propagator 
of a scalar field. 



III. DIRAC OSCILLATOR IN THE (1+1) DIMENSIONAL CASE 

One of the most important and well studied systems in non-relativistic quantum mechanics 
is the harmonic oscillator. This system is characterized by its simplicity and usefulness. The 
harmonic oscillator is a simple system due to the fact that its Hamiltonian operator is quadratic 
both in coordinates and momenta. An analogous system in relativistic quantum mechanics is 
given by a linear interaction term introduced in the Dirac equation. In fact the Dirac equation 
corresponds to the linearization of relativistic Schrodinger equation. So the harmonic potential 
must be introduced in the Dirac equation as the quadratic root of the quadratic potential, i. e. 
the potential must be linear. The interacting system constituted by a relativistic massive fermion 



under the action of a linear potential is known as the Dirac oscillator 

In this section we will consider the Dirac oscillator in the (1+1) dimensional case. Initially 
we will obtain the energy spectrum and the wave functions describing the quantum states of this 
system. Next we develop the canonical quantization procedure for the Dirac oscillator field in 1+1 
dimensions. 



A. Spectrum and wave functions in (1+1) dimensions 

The linear interaction over a relativistic fermion of mass m that moves in one-dimension over 
the z axis is introduced into the Dirac equation 0,0, Q by substituting the momentum p z as 
Pz — > pz — imu^z, where a? is the frequency associated to the oscillator. Thus, the one-dimensional 
Dirac oscillator equation takes the form [h = c = 1) 

d 

The solution of equation allows to obtain both the energy spectrum and the wave functions 
describing the quantum states of this system. As we will show furtherly, these wave functions will 
be used in the canonical quantization procedure of this system. Given that the interaction does 
not mix positive and negative energies, it can be possible to rewrite the four-component spinor 
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state in eq. (jH) as 

/U> ) , (5) 

Jx>/ 

where \4>) and |%) are spinors. If we use the standard representation of the Dirac matrices, the 
time independent Hamiltonian equation takes the form 



m a 3 -(p z +imuz)\ j \<f>)\ ^ I \<j>) 



^3 ■ (Pz - imuz) -m J \\x) J \\x) 

which leads to the following coupled equations 

(E - m) \4>) = 03 • (p z + imujz) \x) , (7a) 
(E + m) \x) = cr 3 • (p z - imujz) \(f>) . (7b) 

Starting from (|7a|) and ([7b]), we obtain that 

[pl + (mu}) 2 z 2 ) |^) = {{E 2 - m 2 )l + mu$) |V) , (8) 

where the commutation relations between the operators p z and z have been taken into account. In 
the last equation, I is the identity matrix 4x4 and j3 is the Dirac matrix 7 . Using the following 
reparametrization ?j] 

C, = (muj)2z, (9a) 
p c = (mw)"2p 2 , (9b) 

it is possible to write the last equation as 

(^ + C 2 )l<A> = ^>, (10) 

where 

E 2 -m 2 

V± = ±1- 11 

raw 

We note that the equations (|7ah and (|7b|) . for the states \4>) and allow to the equation (jlOD 
that has the form of a harmonic oscillator equation. This result suggests us that we can introduce 
the creation and annihilation operators given by 

at = -L(C-^ c ), a=-L(c^ C ), (12) 
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with the purpose to obtain the eigenvalues of the system. In this way, using these operators, the 
equations (|7a|) and (|7b|) have the form 



. y/2moj 

1 (Ton 1 

E — m 



£7 3 a T \x) , (13a) 



E + m 

Substituting again ()13ap into (j!3bp . it is possible to find the following equation for the state \4>) 

. ,. 2mw 



while the state \x) satisfies 



E 2 — m 2 



2muj 



N\<l>), (14) 



Ix> = -=a -i(N + l)\ X ), (15) 

where N = aJa is the occupation number operator of the state. Since N determines the energy 
level of the state where the particle (antiparticle) is, then the energy spectrum can be deduced 
from (1141) and (1151). Therefore we have 



• Prom (|13a|) . the energy spectrum for positive energy states is 

El = 2nmuj + m 2 . (16) 

• From (|13bj) . the energy spectrum for negative energy states is 

E%, +1 = 2(n + l)mu + m 2 . (17) 

These energy spectrums for fermions and antifermions can be written simultaneously if we impose 
the following quantum number condition 

n' = n - 1, with n / 0, (18) 

thus the total spectrum can be written as 

E n = ± {2\n\muo + m 2 Y , con n G Z. (19) 



The upper sign in this expression is taken for n > 0, while the lower sign is for n < |15I |. 
therefore E_ n = —E n (for i. e. the negative quantum numbers correspond to the negative 

energy states. We observe from (I19p . that the lower positive energy state whose energy value is m 
corresponds to the state with n = 0, while the greater negative energy state whose energy value is 



— (2mtjj + m 2 ) 5 corresponds to the state n = — 1. Additionally, we observe that if uj <C m, then 
the energy difference between the states 4>q and X—i is = 2m+w. For a? = 0, i. e. the harmonic 
potential vanishes into eq. then AE = 2m, which is a well known result obtained from the 
Dirac equation in the free case. 

By using the previous results, the states of the system can be written as 



IV>n> = , (20) 




where the quantum number n, which is an integer number, can describe positive and negative 
energy states. Using the expression (Il3b|) . we can write that 



E+m 



a 3 a 



If we apply the occupation number operator N over the state of a system described by \ip n ), we 
obtain 

N\^n)=( I, (22) 

\(\n\-l)\Xn)J 

where we have assumed that the state \4> n ) has an occupation number given by \n\ and where we 
have used the relation (Il3b|) and the properties of the creation and annihilation operators. For the 
last expression, we realize that the lowest spinor \xn) has associated the occupation number given 
by \n\ — 1. Thus the states \4> n ) and \xn) can be written as 

\<j, n ) = \n)tl, (23a) 
|x„> = |n-l>£n, (23b) 

where £n 2 represent the two-component spinors and \n) represents a state with occupation number 
|n|. In consequence, the states of the system are rewritten as 

hM=( (24) 

\\n-i)e n ) 

Taking into account that the creation a) and annihilation a operators satisfy that a) \n) = 
\J\n\ + 1 \n + 1), a \n) = \f\n\ \n — 1), then these operators acting on the state \ip n ) allow to 
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,^TT|n + l>^ +1 



« T l^n) = | V _ / " I , for n / -1, (25a) 

1+1 



v ^ ~ ^ 

whereas these operators acting on the states |^o) y |Y'-l)> which have associate the occupation 
numbers n = 0, — 1, respectively, allow to 

a = \ (1 - P) , (26a) 
a t |^-i) = V^|^o>. (26b) 

In figure Q] we have schematically represented the action of the creation and annihilation operators 
on the positive and negative energy states. We observe that the effect of the annihilation operator 
on the state |V>o) ; which corresponds to the lowest positive energy state, is such that it does not 
annihilate the state but it drives it to the state \tp-i), which corresponds to the greater negative 
energy state. Likewise, the effect of applying the creation operator on the state is such that 

it does not annihilate that state, but drives it to the state l^o)- Therefore, we observe that the 
appearing of the negative energy states generates the well known problem of the Dirac theory: a 
minimal energy state does not exist, then it is possible to obtain an infinite energy amount from 
this system. In order to give a solution to this problem, it is necessary to introduce the Dirac's sea 
picture for the Dirac oscillator which will be performed by means of the canonical quantization 
for this system. 

After calculating the energy spectrum of the one-dimensional Dirac oscillator, we proceed to 
obtain the wave functions. We substitute ([5]) into (|10p . then we obtain the following differential 
equation for the wave function associated to the bispinor \<fi) 



2 



HO = 0, (27) 



where we have used the coordinate representation of the wave function given by 4>(C) = (CI \ 4>)- 
The differential equation (1271) corresponds to the one of a relativistic harmonic oscillator, whose 

n 

solution is [15J 

M0 = N lnl H lnl (0e-^e n - (28) 
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E 



E3 
E 



a \ip 



E-x 

E-2 
E-3 



Figure 1: Energy spectrum of the one-dimensional Dirac oscillator and the action of the creation and 
annihilation operator on some system states. 



Likewise the solution to the differential equation associated to the wave function Xn(C) = (CI \Xn) 
has the form 



Xn(C) = W W -l#|n|-l(C)e~£& 



(29) 



In the expressions for 4> n (C) an d Xn(C)> H n (Q represent the Hermite polynomials. Now, from the 
expression (|13b|) we can obtain the following relation between the spinors and £^ 



e 2 

Sn 



lE n -m x 
-m— cr 3 £ n , 



(30) 



E n + m 

where we have used the properties of the creation and annihilation operators and the definitions 
given by (f23|) . Now it is possible to write 



e n = r 2En , (31a) 




(31b) 



We observe that the spinor £^ is annihilated for the case n = 0, which implies \xo) = 0- So the 



most general solution for the one-dimensional Dirac oscillator equation (|4|) is given by 



1a 



(32) 
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where the normalization for the wave functions has been performed and n = 0, ±1, ±2, .... Finally, 
we can write the one-dimensional Dirac oscillator equation Q in its explicit covariant form [?], [3| 

{vfdp -m + a^F^ = 0, (33) 

where a^ v = §[7^,7^] and the field-strength tensor F^ u associated to the harmonic interaction 
presents in this system is given by 

= d^A„ - d v A^ (34) 
where we have defined the fourpotential associated to the interaction as 

= - (2(u ■ x)x IJL - x 2 u fJi ) . (35) 



In the last expression = (mw, 0) is a fourvector that depends on the reference frame. In the 
references [?], Q|, two different physical pictures for the interaction of the Dirac oscillator have 
been considered. 



B. Canonical quantization procedure in (1+1) dimensions 



By means of the previously defined quantities, we introduce the Lagrangian density of the 
system as 



(36) 



which, by using the Euler-Lagrange motion equation, gives the one-dimensional Dirac oscillator 
equation ()33|) . As the Lagrangian density has an explicit dependence on the position coordinate 



x, it is not possible to determine the energy-momentum tensor of the system 
be found that the Hamiltonian density is 

% = ip^(—ia3 ■ (d z + mujj3z) + j3m)ijj. 



331 ] . However, it can 



(37) 



To follow the standard canonical quantization procedure 
commutation relations of the Jordan- Wigner type for the fermion fields 

\i>a(z,t),ip^(z , ,t)^ = 5 a p5(z - z'), 

Up a (z, t), Ipp {z', t) I = 0, 



311 ]. we initially impose the following 



(38) 
(39) 
(40) 
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Using these relations, we obtain the Hamiltonian operator for the Dirac oscillator field operator ip 
in the form 

H = I dz ft(-ia 3 ■ (d z + muj(3z) + /3m)if). (41) 



Considering that the Heisenberg equation for the field operator ip is given by ip(z,t) = 
—i[ijj(z,t),H], then we can obtain that the motion equation for the Dirac oscillator field oper- 
ator is 

iip(z, t) = (-ia 3 • {d z + mojf3z) + (3m) ^(z, t). (42) 

Now we write the Dirac oscillator field operator using the wave functions of the Dirac oscillator (|32j) 
as the base of the expansion. These wave functions are written in terms of the Hermite polynomials 
which represent a complete set of orthonormal polynomials. Thus, the Fourier serie expansion for 
the Dirac oscillator field operator can be written as 

iE n t 



ip(z,t) = K4) n (z,t)e il 

n=— oo 
00 oo 



n=0 n=l 

= 4+(z,t)+$-(z,t), (43) 

where the positive and negative energy contributions have been separated ip±(z,t) and the spinors 
u n {z) and v n {z) have been defined as 

u n {z) = , v n (z) = . (44) 

\xn(z)ej \X-n(z)e, 



If we make use of the anticommutation relations for the fields ([38]), <\39h and (I40p . we can verify 
that creation bn and annihilation b n operators of positive (for n > 0) and negative (for n < 0) 
energy particles satisfy the following anticommutation relations 

{Win} = S n,m, (45) 

b n ,b m }=0, (46) 
[blP m } = 0. (47) 
Using these anticommutation relations, we can find that the Hamiltonian operator is now given by 

oo 

H= E rMn 



n=— oo 
oo 



: E n bib n - Y KbtJ-n, (48) 



n=0 n=l 
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where we have also used the expansion of the field operator (|43p and the properties of the Hermite 
polynomials. We can observe that the eigenvalues of this Hamiltonian operator can take negative 
values without restriction at all, because it is possible the creation of negative energy particles. 
To solve this problem we take the Dirac's sea picture for this system, so we impose that all the 
negative energy states are occupied by negative energy particles. This configuration is called the 
vacuum state of the Dirac oscillator field |0) that is written as 34| 



|0> = IP-|0d>, (49) 

n=l 

where \0d) is the Dirac vacuum state which is characterized for not been occupied by fermions of 
positive or negative energy. We observe that if we apply the creation operator of negative energy 
particle b[_ m |0) on the vacuum state |0) we obtain 

oo 

bl m \0) = l{bl m bl n \0 D }=0, (50) 

n=l 

which implies that it is not possible to create a new negative energy fermion because all the negative 
energy states are occupied. On the other hand, it is possible that a negative energy state can be 
annihilated and then a hole can be originated in the Dirac's sea. This hole represents an antiparticle 
with negative energy. If we use the anticommutation relations (|45|) . (|46p and (|47p into (|48|) . we 
find that the Hamiltonian operator takes the form 



H = J2 E J> n b n + Y.E n (L n 6 f _ n - l) . (51) 



n=0 n=l 

If the divergent energy of the vacuum Xm^=i En is subtracted, the Hamiltonian operator can be 
reduced to 

oo oo 

H = EJlK + E n b- n bl n . (52) 

n=0 n=l 

We note that now this Hamiltonian operator does not have the problem observed before because it 
has been positively defined. Finally, we perform the following canonical transformations to clearly 
differentiate between the operators associated to particles and antiparticles [31 ] 

Hi = 4> ( 53a ) 

b n = dn, (53b) 

S_ n = d\, (53c) 

bL n = d n , (53d) 
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where it is possible to identify d n and d n as the annihilation and the creation operators of an- 
tiparticles, respectively. These operators satisfy the anticommutation relation |d n ,dJn| = 5 n>m . 
These transformations imply that now the application of the annihilation operator of negative en- 
ergy particle on the vacuum state can be understood as the application of a creation operator of 
antiparticle. Using this notation, the Fourier serie expansion for the Dirac oscillator field (143p has 
now the form 

oo oo 

#(z, t) = J2 c n n n (z)e- lE " t + $ n v n {z)e iE ^. (54) 

n=0 n=l 

Now the Hamiltonian operator for the Dirac oscillator field can be written as 

oo oo oo oo 

H = Y, E nc ] n c n + ^ E nd\dn = £ E n n$ + £ , (55) 

n=0 n=l n=0 n=l 

where ri$ represents the number operator for fermions and n$ represents the number operator 
for antifermions. It is possible to observe that cli acting on |0) creates a fermion having energy E n 
and charge +e, whereas d\ creates antifermion having energy E n and charge — e. In this scheme, 
an interacting relativistic massive fermion is described as an excitation from a degree of freedom 
of the Dirac oscillator quantum field. For the Dirac oscillator field the energy quanta E n are the 
energies of linear harmonic oscillators. We observe that, in this sense, the Hamiltonian operator 
(|55p is analogous to the Hamiltonian operator for the free Dirac field in three dimensions given by 

Using the Fourier expansion (I54p it is possible to determine other relevant physical quantities. 



For instance, the charge operator Q = ef dz^ip 31f] can be written as 



Q = e^2 (Hfin ~ 4cn) , (56) 
n=0 

where the non-observable charge of the vacuum has been removed. Likewise, the momentum 



operator P = —i f dzip^Vip 3l|] can be also written in an explicit form using the expansion of the 
Dirac oscillator field. 



The Feynman propagator for the Dirac oscillator field in the coordinate space is defined as 



31 



iS F p {z-z',t-t') = (0\f(j> a (z,t)^(z',t')) |0). (57) 
Substituting (143 p into (|57p and taking into account the definition of the time-ordered operator T, 



15 



we obtain 



iS^(z -z',t- t') = Q(t - t') J2 u a , n {z)u 0>n {z')e- iE ^-^ 

oo 

- Q(t' -t)J2 VoA^pAz'Y^^ ■ 



(58) 



Now we will obtain the Feynman propagator in the momenta space. To do it, we consider firstly 
that the Fourier transformation for the Hermite polynomials can be written as 

(59) 



$U~H n {x)\ = (-i) n e~-H n (k), 
and secondly, that the contour integral for the energy eigenvalues is written as 



dp e tPo ^ * ) 
C 2vr p 2 - p 2 n 



e(t - 1') 



-iE n {t-t') 



2E„ 



+ &(t'-ty 



2 E n 



where p\ = 2 |n| muj + m 2 and E n is given by (|19p . Then we find that 



(60) 



(61) 



c (2^) 3 

where the Feynman propagator in the momenta space for the one-dimensional Dirac oscillator field 
is given by 



Sap{Po,Pz,Pz') = ^2 — 2 [ u a,n{Pz)up,n{Pz') ~ U Q ,n-lfc)^,n-lfe')] 

n=0 P 



(62) 



211 ] . where the Feynman propagator of the one- 



This result is in agreement with the presented in 
dimensional Dirac oscillator was obtained using the path integral formalism and working with a 
different representation of the Dirac matrices. 



IV. DIRAC OSCILLATOR IN THE (3+1) DIMENSIONAL CASE 

In this section we will consider the Dirac oscillator in (3+1) dimensions. Initially we will obtain 
the energy spectrum and the wave functions describing the quantum states of this system. Next we 
develop the canonical quantization procedure for the Dirac oscillator field in the 3+1 dimensional 
case. 



A. Spectrum and wave functions in (3+1) dimensions 

In order to obtain the Dirac oscillator equation in the (3+1) dimensional case, we perform the 



following substitution which is analogous to the minimal substitution 



,G,Q 



p — > p — imuj r Y°r , (63) 
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where p = (p x ,Py,Pz) and f = (x,y,z). The Dirac oscillator equation has the form 

.8 , „ 



'at 



fa • (p — imoj^°r) + /3m) 



(64) 



The solutions of the Dirac oscillator equations are described by the spinor states \ip), that can be 
written as 



(65) 




where \(p) and \\) are bispinors. The bispinors obey the following equations 

(p 2 + (mw) 2 r 2 ) \tp) = (E 2 - m 2 + (3 + 4s • t)mu) \tp) , 
(p 2 + (mufr 2 ) \x) = (E 2 - m 2 - (3 + 4s • L)mw) |x) • 



(66a) 
(66b) 



From these equations, it is evident that the Dirac oscillator in the (3+1) dimensional case presents 

nn 

a strong spin-orbit coupling term [3, [14J . It has been demonstrated that the Dirac oscillator is a 
system where the angular momentum is conserved , thus the quantum numbers of full angular 

momentum j and parity are good quantum numbers. Therefore, it is convenient to separate the 
energy eigenfunctions in two parts: radial and angular. The spatial coordinate representation of 
the states (I65p allows to the spinor wave functions given by 



r \iG n , K (r)y- K , 9 (0i< 



(67) 



where n is the principal quantum number, k is a quantum number related with the angular mo- 
mentum and parity, and g is a quantum number related to the projection of the angular momentum 



in the z axis. The quantum number k is defined as 



35] 



-(Z + l), if j = l + \, 
I, if j = l-\. 



(68) 



The full angular momentum j takes values j = I ± i, where I is the angular momentum and \ is 
the spin of the fermion. The angular momentum V associated to the upper and lower components 
is 



35] 



V = 2j-l 



l + l, if j = l + \, 



(69) 



1-1, if j = l 



l 

2' 
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The spinorial spherical harmonics y K ,g(6,4>) an d 



are given by 35] 



y-n, 9 ie,4,) 



2j+2 2 ± 



j±5±l y 
2J+2 J Z',S+i 




(70a) 



(70b) 



) Qi 



Substituting ((671) into the Dirac equation (j64l) represented in coordinate space, we can find the 

(71a) 
(71b) 



+ 



F n Jr) = (E + m)G n Jr) 



dr 

Tr + J G n , K (r) = {E— m)F n , K (r). 



The solutions of this equation system are BQ.Q 



, j_ii mujr^ /+— , Ov 

Fn,K\T) = 4[vmwr] t+ e 2 L n 2 (mujr z ), 



(72a) 
(72b) 



G n ,K,{r) = ±sgn(K)A'[^mujr] l ' +1 e ^ 2 (mur 2 ) , 

where L l n (x) are the Laguerre associated polynomials and A, A' are normalization constants given 
by 

Imuj |n|! (E n , K + m) 



" _r(|n| + Z + 3/2)£ niK 

^Jmuj \n'\l (E n ^ K — m) 
T{\n'\ +l' + 3/2)E ntK 

where the quantum number \n'\ takes the following values 



(73a) 
(73b) 



n 



\n\ 



11 



(74) 



1, for k < 
for k > 

The radial functions F(r) and G(r) obtained here have the same radial structure as the one 
associated to the non-relativistic three-dimensional harmonic oscillator. On the other hand, the 
energy spectrum for this case depends on the k value. It is possible obtain that the energy 
eigenvalues are [lj]] 



For k < 



E n ,K = ± 



n\ muj, 



where the quantum number n can take the values n = 0, ±1, ±2, . . 
chosen for n > 0, meanwhile the negative sign is chosen for n < 0. 



(75) 

and the positive sign is 
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For k > 

E n ,K = ± ^rn 2 + 4 ^|n| +7+ mw, (76) 

where n = ±0, ±1,±2,... and the sign is chosen as was mentioned in the before item. 
Moreover, the energy spectrum satisfies the following symmetry condition E- U)K = —E n}K , 



except for n = which means k < 0. For t 
the quantum numbers n = +0 and n = — 



lis case, it is necessary to differentiate between 



1 



With the purpose of implementing the canonical quantization procedure for the Dirac oscillator in 
(3+1) dimension, we observe that the spinorial wave functions (|65p satisfy the following orthonor- 
mality and completeness relations [lj| 

OO 2 oo 

E E E i>n,.A^U 9 (r) = 5 3 (?-r)h, (77a) 

^f^L^gi^n'^^g'ir) = £n,n'<W'<W , (77b) 

respectively. Finally, we note that the covariant Dirac oscillator equation for the (3+1) dimensional 
case has the same form as the equation (|33p for the (1+1) dimensional case. 



B. Canonical quantization in (3+1) dimensions 

To proceed with the canonical quantization for the Dirac oscillator field in the (3+1) dimensional 
case, we observe that the form of the Lagrangian density for this case is similar to the one in the 
(1+1) dimensional case and it is given by 

L = 4>(i^d^ - m)^ + Fp,. (78) 

As in the (1+1) dimensional case, the Lagrangian density depends explicitly on the spatial coor- 
dinates and then it is not possible to obtain the energy-momentum tensor of the system. How- 
ever, by means of a Legendre transformation, we can obtain the following Hamiltonian density 
% = i/j^(—icx ■ (V + mujfir) + j3m)ip. Thus, the Hamiltonian of the field is 

H = d 3 ri[) j (-ia ■ (V + mw^r) + /3m)%/j. (79) 
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To perform the canonical quantization for this case, as is usual, the fields are considered as field 
operators. Then we establish the following Jordan- Wigner commutation relations 

{^(r,i),^,t)} = <^ 3 (r-fO, (80) 
{^ a (r,i),^(r^)}=0, (81) 
{&(r,*),$(»V)}=0, (82) 

that allow us to write the Hamiltonian operator as 

H= j d 3 r$(-ia ■ (V + mujf3r) + Pm)$. (83) 

Now we can expand the field operator ip(f,t), using the spinor wave functions (|65p . as follows 

OO 2 oo 

^*)= E E E ^>»«( r >" ,£ "'" (84) 

where the operators b n ^,g and Sn )Kj g, respectively, annihilates and creates fermions in a state defined 
by the quantum numbers n, k and g. These operators obey the following anticommutation relations 

{Vre,-?,^^/} = S n,n' S K,K'Sg,g>, (85) 

y>n,K,g, b n i K i g if = 0, (86) 

Using the properties of the energy eigenfunctions (|77ap and ()77b|) into the Hamiltonian operator 
(|84p . we can obtain 

oo 5 oo 

^ = E E E E n,Mi,K,gK,K,g- (88) 
k=— oo „— n=— oo 

This operator can be rewritten by splitting the positive and negative energy contributions and 
taking into account the two types of spectrum depending on the k value. In this way, we obtain 



oo oo oo oo 



H — ^2 ^2 ^2 E nt J}\ l K g b n ^^ ^2 ^2 ^2 E n ,Kb- n ,K,gb-n,K,,g 
k=1 g n=0 k=1 g n=0 

oo oo oo oo 

+ E E E E n-nb n -K,gbn-K,g ~ E E E E n-3>- n ,- K , g b-n-K,g, (89) 
k=1 g n=0 k=1 g n=l 

where we have used the convention E n _ K = E n ^ K) for k < 0, and we have omitted the limits of the 
sum over g. It is evident that this Hamiltonian is not defined positively as happened in the (1+1) 
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dimensional case. Again we use the picture of the Dirac's sea in order to solve this problem. To 



do it we take the vacuum state in an analogous way as it was defined in 

oo oo 

i°)= n nn s -n,«jM- (90) 

k=— 00 g n=0 

Therefore we can rewrite the Hamiltonian operator without considering the vacuum energy in the 
following way 

00 00 00 00 

H' = E E E E ^iMi,n,gK,K,g + E E E E n,Kb-n,K,gb^ n ^ g 
k=1 g n=0 k=1 g n=0 

oo oo 00 00 

+ E E E E n,-X,- K>g k,-K, 9 E E E E n,-j>-n,-K,gbl n ^ g . (91) 
k=1 g n=0 k=1 g n=l 

We observe that this Hamiltonian operator describes two different types of particles and antipar- 
ticles, because particles and antiparticles have different energy spectrums depending on the sign 
of the k value. This value depends explicitly on the full angular momentum j due to the value 
of j is based on the spin state of the fermion. In this way, we perform the following canonical 
transformations 

(92a) 
(92b) 
(92c) 
(92d) 

where 6n,K,g is the creation operator of particles with full angular momentum j = I — ^, i. e. 





K.g u n 


k.,9 




K,g — °n 


K,g 


K- 


K,g ™"n 


,K,g 


-n- 


n,g = fn 


,K,g 



particles with k > 



328; Cn,K,g is the creation operator of antiparticles with full angular momentum 



j = I — i, i. e. antiparticles with k > 



32]; d n .K,,g is the creation operator of particles with full 



32]; fn, K ,g is the creation operator of 



angular momentum j = / + ^, i. e. particles with k < 

antiparticles with full angular momentum j = I + ^, i. e. antiparticles with k < 

Finally, we can obtain that the field operator for this case can be written by means of the 
expansion 

oo oo oo oo 

*) = E E E V^n,«>>-^ + E E E 4,^-n, K >>^< 
k=1 g n=0 k=1 g n=0 

oo oo oo oo 

+ E E E d n ,^n,- K ,g(^- iEn -- Kt + E E E fL^-n,- K>9 (?)e iEn '- Kt - (93) 

k=1 g n=0 k=1 g n=l 

In an analogous way as was showed for the (1+1) dimensional case, starting from the expansion for 
the field operator (|93p . it is possible to obtain the different relevant physical quantities associated 
to the Dirac oscillator field in (3+1) dimensions. 
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V. CONCLUSIONS 



In this work we have performed the canonical quantization of the Dirac oscillator field in (1+1) 
and (3+1) dimensions. This quantization has been possible because the solutions of the Dirac 



111 ]. If this fact would not have satisfy, then 



331 ]. which had 



oscillator equation do not present the Klein paradox 

it had not been possible to distinguish between positive and negative energy states 
restricted the possibility to perform a Fourier expansion of the field operator. The Dirac oscillator 
field has been quantized following a similar procedure as if this field were free 3l[. However, 
this procedure implies differences with respect to the Dirac free field quantization. For instance, 
the Dirac oscillator is an interacting system in which there exist bound states define with specific 
quantum numbers {{J. In this case the field operators create and annihilate fermions with well 
determined energy values in contrast to what happens in the free field case where the states have a 
well defined momentum [22]. Moreover, for this case, the momentum operator depends explicitly 
on the time thus the fermions created do not have a determined momentum. We have found that 
for the Dirac oscillator field it has been impossible to obtain the energy-momentum tensor due 
to the fact that the Lagrangian density has an explicit dependence on the spatial coordinates. 
Nevertheless, this problem could be solved by the introduction in the system of an additional field 
which describes the interaction. We have also obtained the Feynman propagator in the (1+1) 
dimensional case in agreement with the result obtained in the literature by using an functional 
procedure 



211 ]. We note that for the Dirac oscillator field in the (3+1) dimensional case, we have 
found that there exist two types of particles and antiparticles because there are two possible values 
for the full angular momentum. Finally, we have found that while for the free Dirac field the energy 
quanta of the infinite harmonic oscillators are relativistic energies of free particles, for the Dirac 
oscillator quantum field the quanta are energies of relativistic linear harmonic oscillators. The 
canonical quantization procedure for the Dirac oscillator field in (2+1) dimensions is an exercise 
to develop explicitly. We consider that the possibility to study the Dirac oscillator as a quantum 
field opens the doors to future applications in different areas of the Physics. 
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